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The glueball spectrum from novel improved actions. 
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Results for the inter-quark potential and low- lying SU(3) glueball spectrum from simulations using a new 
improved action are presented. The action, suitable for highly anisotropic lattices, contains a two-plaquette term 
coupling with a negative coefficient as well as incorporating Symanzik improvement. 



1. INTRODUCTION 

The QCD glueball spectrum has been inves- 
tigated in low-cost simulations using anisotropic 
lattices 0. To reduce the computational over- 
head, the spatial lattice was kept rather coarse 
(0.2-0.4 fm) while the temporal spacing was made 
much finer. The fine temporal grid allows ade- 
quate resolution of the Euclidean-time decay of 
appropriate correlation functions which, for glu- 
onic states are rather noisy and fall too rapidly 
on coarse lattices. 

In these simulations, the scalar glueball suf- 
fered from large finite cut-off effects. The mass 
in units of tq fell sharply until the spatial lattice 
spacing, a s was about 0.25 fm when the mass rose 
again; the "scalar dip" . 

At the conference last year, we presented re- 
sults from simulations with an anisotropic Wil- 
son "two-plaquette" action which included a term 
constructed from the product of two parallel pla- 
quettes on adjacent time-slices ||. This was 
found to reduce the scalar dip significantly. Here, 
we report on the status of simulations in progress 
using a Symanzik-improved action including a 
similar two-plaquette term. 

In this study, we tune the anisotropy parame- 
ter in the lattice action to recover Euclidean in- 
variance in the "sideways" potential. With these 
parameters fixed, we investigate the inter-quark 
potential for this action as an initial test that the 
benefits of the Symanzik program are preserved 
by the addition of the extra term. We are cur- 



rently computing the glueball spectrum for this 
action. 

2. THE ACTION 

Following Ref. |], we begin with the plaquette 
operator, 

Pyvix) = ^ ReTr U»(x)U v {x+ii)Ul(x+v)Ul(x). 

The Wilson (unimproved) discretisation of the 
magnetic field strength is then constructed from 
the spatial plaquette. 
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where i,j are spatial indices and £o is the 
anisotropy, a s /a t at tree-level in perturbation 
theory. We introduce a term which correlates 
pairs of spatial plaquettes separated by one site 
temporally 

^s t] = \ E {l-^i(z)^i (* + *)}• (2) 

The separation of the two plaquettes allows the 
standard Cabibbo-Marinari and over-relaxation 
gauge-field update methods to be applied. In- 
cluding two-plaquette terms adds a computa- 
tional overhead of only 10% to our improved ac- 
tion workstation codes. 

It can be shown that for all ui, the operator 
combination, 
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Figure 1. The ratio of Eqn. || for n = 3, 4, 5 



has an identical expansion in powers of a S) t (at 
tree-level) to f2 s up to 0(af). Thus, starting from 
the improved action Sn used in Refs. Jl],[|, it 
is straightforward to construct a Symanzik im- 
proved, two-plaquette action by simply replacing 
the spatial plaquette term in Sn with the linear 
combination tl s of Eqn. |[ In full, this action is 



I 



5(l+a;) 
3uj 
4 

3*? 



fi.. 



fit 



5w 
3^f 



fii 2t ) 



12u° 



rfi^ 



1 



fij 



(4) 



with f2 t the temporal plaquette and fli R \fl[ R ^ the 
2x1 rectangle in the (i, j) and (i, i) planes respec- 
tively. This action has leading 0(af , a 2 , a s a 2 ) 
discretisation errors and only connects sites on 
adjacent time-slices, ensuring the free gluon prop- 
agator has only one real mode. 

The free parameter w is chosen such that the 
approach to the QCD continuum is made on a 
trajectory far away from the critical point in the 
plane of fundamental-adjoint couplings. Close to 
the QCD fixed point, physical quantities should 
be weakly dependent on u). This provides us with 
a consistency check, however the data presented 
here are for one value only, u) = 3. 
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Figure 2. The inter-quark potential of Sn with 
a s « 0.22 fm for sources separated along different 
lattice axis. 



3. TUNING THE ANISOTROPY 

At finite coupling, the anisotropy measured us- 
ing a physical probe differs from the parameter in 
the action at 0(a s ) ||- In previous calculations, 
we relied upon the smallness of these renormal- 
isations for the (plaquette mean-link improved) 
action Sn . For the action of Eqn. |[ these renor- 
malisations are larger and thus we chose to tune 
the input parameter in the action to ensure that 
the potentials measured along anisotropic axes 
matched. We follow a similar procedure to Ref. 
||. The potentials between two static sources 
propagating along the z-axis for separations on 
both fine and coarse axes, V t and V s respectively, 
are measured using smeared Wilson loops. Since 
the UV divergences due to the static sources are 
the same, tuning £o such that the ratio 



a s V s (na s ) 
a s V t (mna t ) 



1. 



(5) 



implies the anisotropy £y = m (m G Z). A con- 
sistency check is provided by studying different 
coarse source separations, na s . Fig. [I] shows 
this tuning for n = 3,4,5, where the desired 
anisotropy is 6. Consistency is observed for n = 4 
and 5 and the appropriate £o is found to better 
than 1%. 
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Fi gure 3. The scslar glueball of both Symanzik Figure 4. The scalar, tensor and pseudoscalar 

improved actions, S u and S u glueballs of 5 W 



4. SIMULATION RESULTS 

4.1. The inter-quark potential 

The replacement of the spatial plaquette in 
Sn with fl s of Eqn. |^ should lead only to 
changes in the irrelevent operators responsible for 
0(a A s ,a s a 2 s ) errors. To test this replacement still 
generates an improved action with the good ro- 
tational invariance of Sn , the inter-quark poten- 
tial was computed for a variety of different inter- 
quark lattice orientations. The potential is shown 
in ||, and shows excellent rotational invariance. 
We conclude that the benefits of the Symanzik 
improvement programme are preserved by includ- 
ing the two-plaquette term for a typical value of 
oj useful for glueball simulation. 

4.2. The glueball spectrum 

At present, we are computing the glueball spec- 
trum on the £y = 6 tuned lattices. Preliminary 
data are presented in Figs. || and [|. In Fig. 
|[ the finite-lattice-spacing artefacts in the scalar 
glueball mass for the new action are compared 
to those of Sn . The lattice cut-off dependence is 
seen to be significantly reduced and for the range 
of lattice spacings studied here, the mass rises 
monotonically with lattice spacing rather than 
falling first to a minimum. Fig. [| shows the lat- 
tice spacing dependence on the tensor and pseu- 



doscalar glueballs. Their lattice spacing depen- 
dence is similar to the form for Sn and consistent 
with leading 0(a A s ) behaviour. 

5. CONCLUSIONS 

Preliminary data from our simulations of the 
Symanzik improved action of Eqn. || suggest 
the scalar dip is removed by inclusion of a two- 
plaquette term with negative coefficient, consis- 
tent with the argument that the poor scaling of 
the scalar glueball, even after Symanzik improve- 
ment, is caused by the presence of a nearby crit- 
ical point. 

The inter-quark potential on this new action 
exhibits equally good rotational symmetry to the 
improved actions of Refs. 
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